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Abstract 

We investigate the relationship between the effective diffusivity and effective 
drift of a particle moving in a random medium. The velocity of the particle 
combines a white noise diffusion process with a local drift term that depends 
linearly on the gradient of a gaussian random field with homogeneous statistics. 
The theoretical analysis is confirmed by numerical simulation. 

For the purely isotropic case the simulation, which measures the effective drift 
directly in a constant gradient background field, confirms the result previously 
obtained theoretically, that the effective diffusivity and effective drift are renor- 
malized by the same factor from their local values. For this isotropic case we 
provide an intuitive explanation, based on a spatial average of local drift, for the 
renormalization of the effective drift parameter relative to its local value. 

We also investigate situations in which the isotropy is broken by the tenso- 
rial relationship of the local drift to the gradient of the random field. We find 
that the numerical simulation confirms a relatively simple renormalization group 
calculation for the effective diffusivity and drift tensors. 



DAMTP-96-34 



1 



1 Introduction 



A problem of great interest concerns the computation of effective parameters for a 
diffusion process that combines molecular diffusion with a drift term linearly depen- 
dent on the gradient of a random scalar field with spatially homogeneous statistics. In 
the case where the system, including the statistical properties of the random field, is 
isotropic it has been shown, on the basis of plausible but not completely established 
assumptions, that the effective long range diffusivity is accurately predicted by a renor- 
malization group calculation (RGC) |j]> @> SI- The resu lt) which is embodied in a very 
simple formula, is confirmed by numerical simulation to high accuracy. The result is 
also consistent with the straightforward perturbation expansion to two-loop order but 
remains accurate beyond the applicability of the expansion to this order |fj. |4| . 
An important parameter relevant to the long time behaviour of the system, that has 
not generally been studied, is the value of the effective drift parameter, or tensor, that 
determines the average drift of particles in an imposed large-scale field gradient. A fea- 
ture of the isotropic case, respected by the RGC, is that this effective drift coefficient 
is renormalized relative to its local value by the same factor as the effective diffusivity 
is renormalized relative to the molecular diffusivity |T|, [| |J. By including an ap- 
propriate constant drift term in the original molecular diffusion process we confirm, in 
this paper, directly by numerical simulation, the equality of the two renormalization 
factors. We also give an intuitive explanation in terms of spatial averages, of why one 
should expect that the drift coefficient be renormalized. 

Isotropy can break down for three reasons. First, the statistics of the random field may 
not be isotropic. We have investigated this in a previous paper || and shown that, 
although the predictions of the RGC are still reasonably accurate there are emerging 
disparities with the results of the numerical simulation. Such an outcome is to be 
expected since the results of the RGC were, in this case, shown to differ from straight- 
forward perturbation theory at two-loop order. Second, the molecular diffusion tensor 
may be non-isotropic and third, the same may be true of the local drift tensor. It is this 
third possibility that we investigate in this paper by numerical simulation and compare 
with the predictions of the appropriate RGC. The results actually agree rather well 
suggesting that the isotropy of the random field statistics is an important condition for 
the success of the RGC. 
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2 Diffusion with Background Drift 

It is convenient for the purpose of discussing anisotropy to consider the most general 
diffusion equation of the class in which we are interested. It has the form 

|| = ft«a,--A9.^(x)-u i )P . (1) 
Here P is the probability density of a particle moving according to the equation 

Xi = A°.^-0(x) +u t + w t {t) , (2) 
where the white noise term has zero mean and the correlation function 

(w l {t)w,{t')) = 2K%5{t-t l ) , (3) 

and ui is a constant drift term. In appropriate units we would expect the drift to be 
given by 

Ui = ^j93 > ( 4 ) 

where gi is a uniform gradient field. 

The equation for the static Green's function corresponding to a unit source at x' is 

d t (k%8 3 - A?.fi^(x) - «i) G(x, x') = -*(x - x') . (5) 

The effective Green's function obtained after averaging over the random ensemble of 
flows we denote by C?(x — x') . It is simpler to study this function in terms of its Fourier 
transform. It has the form 

Q(k) = K n k m k n - S(k) + m^-(k)]- 1 . (6) 

At small k the irreducible two-point function S(k) satifies 

E(k) ~ Gijkikj , (7) 

with the result that the effective long range diffusivity is 

K ij = K ij ~ a ij ) (8) 

and 

Wj(k) ~ kifXij. (9) 



for some coefficient //^ that we now evaluate. We show below that 

W i — [(^°)~ 1 ]imKn(0, k) , (fO) 

where V^(q, k) is the (Fourier transform of the) vertex function that measures the 
influence of a weak external field on £(x, x') . It is calculated below. For small k we 
have 

V m {0, k) ~ k n X nm , (II) 

where \ mn is the effective coupling referred to in the introduction. It follows that for 
small k 

(12) 

The interpretation of this result is that the effective drift is 

«™ = A m „5f n . (13) 

The measurement of the effective drift in a given external field allows us to extract the 

effective coupling Ay from the simulation. 

For the purposes of simulation we assumed that 

a(x-x') = / ^^(q)^ q - (x - x,) , (14) 

with 

D(q) = ^e^W . (15) 
k 

With this normalization 

<(0(x)) 2 > = A(0) = 1 . (16) 

3 Graphical Rules for Perturbation Theory 

The Feynman rules for the diagrammatic perturbation expansion are essentially tha 
same as in the isotropic case. We have 

(i) The sum of the inwardly flowing wave-vectors at each vertex is zero. 

(ii) Each solid line carries a factor of 1 , • 

G^q 

(iii) Each loop wave- vector q is integrated with a factor . 



(iv) Each vertex of the form of Fig. 1 carries a factor ftA^- (k + q) .. 

(v) Each vertex of the form of Fig. 2 carries a factor —iujkj. 

(vi) Each dashed line carries a factor D (q). 

4 One-Loop Contributions 

The one-loop contribution to S(k) is associated with Fig. 3 . It is, according to the 
Feynman rules, 

L (k) " J (2nr D[q) (k + qK s (k + q) s ' { ' 

The one- loop contribution to WjWi(k) is associated with Fig. 4 . The Feynman rules 
imply that it has the form 

H/ / n [ d 3 q (k + q)^ jqj (-iUjik + qjj) k m \° mn q n 
UiWAK.) = — / r~7r — r : — T—r, r^r; ■ (18) 

If we compare this result with that for the general vertex at one-loop associated with 
Fig. 5, namely 

V-( k') - f rf3q ^ + gji^g A ^ k ' + ^ ^rrXmQn , , 

ilq ' } ~ J (27r)3(k + q) r K0 s (k+q) s (k' + q) u <(k' + q), ' 1 J 

we see immediately that 

^(0,k) = A°^(k) . (20) 

This result is easily generalized to all orders in perturbation theory. It is equivalent to 
the relation used in a previous paper f|. 

5 Spatial Averages and Mean Drift 

Some insight into the renormalization of the drift coefficient in the isotropic case (A^ = 
XoSij and n\f = KoSij) can be obtained by considering a situation in which we allow a 
large cloud of particles to equilibriate in a large cubical sample with periodic boundary 
conditions. We then disturb the situation by imposing a uniform background drift field 
g causing a mean drift in the particle cloud. The average velocity of the particles is 
then renormalized relative to the mean local drift because of the lack of uniformity in 
the background probability distribution adopted by the particles. 
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In the absence of a background drift field the particles adopt a density distribution 
Po(x) that yields zero particle current. That is 

J = ( Ko V - A o V0(x))P o (x) = . (21) 
The solution of this equation is 

P (x) = iVexp{%(x)} , (22) 
where we will choose N so that Po(x) can be viewed as a probability distribution. 

/cube d3xfWx) = 1 ' (23) 
For a very large cube it is acceptable to enforce this equation as an ensemble average. 
This allows us to evaluate iV as 

ff= ( vop (KS' A(0) )) • (24) 

where V is the volume of the cube. 

When an external constant drift g is applied the system equilibriates with a new 
distribution P + Pi that satisfies 

V.(/toV-AoV0(x)-A o g)(Po(x)+Pi(x)) = . (25) 

It turns out that no change of normalization is required. This leads to the solution, 
correct to 0(g), for Pi 

P 1 (x) = -A |rf 3 x / G'(x,x / )VPo(x / ).g • (26) 

where G(x, x') is the Green's function for the problem without drift. 
At a point x in a given sample of the medium the velocity of a particle is A o (V0(x) + g) 
after averaging over molecular diffusion effects. The the drift u of particles in the steady 
state situation with drift can be obtained as a spatial average of this velocity over the 
large cubical sample. 

u = A J rf 3 x(V0(x) + g)(P (x) + P x (x)) (27) 

To 0(g) we have 

u t = Xotfa -Xoj rf 3 xrf 3 x , (<9^(x))G(x,x')^Po(x')](? J . (28) 



In fact it is reasonable now to take the ensemble average of this quantity, with the 
result 

u t = Xo[5ij -Xoj C / 3 xrf 3 x'((^0(x))G(x,x / )^P o (x'))]^ • (29) 

To lowest (non-trivial) order in Ao we can replace G(x, x') by Go(x — x') and use the 
identity 

(0(x)P o (x')) = -A(x-x')/^ , (30) 

K 



to obtain 

A 2 

u 

Expressing this in Fourier space we have 



Xo[Sij - ^ / rf 3 xG (x)9^.A(x)]^ . (31) 



Taking into account the isotropy of the statistical ensemble that we assume in this 
case, we see that u = A e g where 

Ae = Ao (l - ^ A (°)) ( 33 ) 

This is identical with the one-loop perturbation result of previous work |L|, 0, [3|, |j. 
A careful analysis of the two-loop diagrams confirms the result to this order. This of 
course, is consistent with the RGC result 

Ae = Aoexp|-l^|A(0)| , (34) 

discovered previously. From this point of view then the origin of the renormalization of 
the effective drift term comes about because in a steady state the particles adopt, in a 
given sample, a non-uniform distribution, appropriate to the sample, and the resulting 
spatial average of the local drift is modified relative to the local ensemble average 
which yields the local mean value. This should be contrasted with the situation in 
incompressible flow where the steady state distribution of particles inevitably remains 
uniform leading to a spatial average that coincides with the ensemble average and no 
renormalization of the drift coefficient. 



6 Simulation of Drift in the Isotropic Case 

In Table 1 we exhibit the results of measuring both and Ae, for an isotropic situation, 
over a range of values of the disorder parameter Ao/^o • We assume k = 1 and k = 1 
throughout. The simulation, 
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N 


Ao 


9 


«e 


(Ae/Ao) 


128 


1.0 


.05 


.722(2) 
(.717) 


.726(2) 
(.717) 


256 


1.5 


.05 


.470(5) 
(.472) 


.474(2) 
(.472) 


256 


2.0 


.05 


.269(4) 
(.264) 


.266(11) 
(.264) 


512 


2.5 


.05 


.139(5) 
(.125) 


.140(20) 
(.125) 



Table 1: N is the number of modes in random field. Theoretical values in brackets. 

which is based on a continuum construction for 0(x), has been described in previous 
papers |3|, |], [5]. An important point for the present paper is that the random field is 
constructed from a set of N modes and the integrity of the Gaussian property of the 
statistics of the random field is dependent on having a sufficiently large value of iV . 
The results clearly show the equality of the two renormalization parameters for a wide 
range of disorder parameter. This common value is also equal, as was found in previous 
work, to the RGC prediction of exp { — |Aq/kq}. There is a slight discrepancy at the 
higher values of the disorder parameter. We feel that this is a systematic error in 
the simulation due to the limited number of modes incorporated in the random field. 
Another possible source of error is that the value of the drift parameter has become 
so large that 0(g 2 ) effects are influencing the values of the measured quantities. It is 
also possible that the assumptions behind the renormalization group calculation may 
no longer be valid at these values of Ao • It is interesting to note that, nevertheless the 
equality of the two renormalization factors is maintained throughout with particular 
accuracy. We can be reasonably confident therefore of our renormalization group results 
and the associated Ward identity J|, || in this isotropic case. 

7 "Ward" Identity 

In previous work we suggested a Ward identity as an explanation of the proportionality 
of the renormalization of the vertex and diffusivity matrix. The identity was verified to 
two-loop order in perturbation theory. Here we show that this Ward identity changes 
form when the bare vertex and diffusivity matrices are no longer proportional to one 
another. We will only discuss the the change at one-loop order since this is sufficient 



8 



to demonstrate the breakdown. 

From the formulae of the previous section we see that to one-loop 

A E(k ) = [^-i] ..y. (k; k) + Ui(k) > (35) 

where o o o o 

<( J ~ J (q + k^q + k), • (36) 

For small k we see that 

W - - / |^(D - ^^.^ } + o(^ 4 ) ■ M 

It is easy to see that that the right side of this equation vanishes to 0(A; 4 ) when k° oc A . 
In general however it does not do so. A simple case to consider is one in which the 
statistics of the 0(x)-field are isotropic, the diffusivity has the form «;?. = kqS^ but the 
drift coefficient retains an anisotropic tensorial structure. We easily evaluate the right 
side of eqflTFD to be 

2A(0) [3 ; 
3k I - 

If we define 

1 

3' 

so that A is the traceless or quadrupole part of A then we can recast the above 
equation in the form 

2A(0) f3. 
3k 

that is 



= {^[(A ) 2 ]^- + \\ mm k\%k 3 } + Oik') , (40) 



kiUiOc) = -^^-k^XXjkj + 0(k 4 ) . (41) 

This form of the result makes it clear that when the traceless part of the drift tensor 
A vanishes we return to the situation in which the original Ward identity holds. 
For the particular case we are considering the the modified Ward identity implies for 
small k the result 

w m = ~tB ( a - a -- + 2A ^} kj ~i^r {l [m ^ - \ x ^%k 3 }+o{k*) . 

(42) 

On the basis of this calculation we do not expect a simple relationship between the 
macroscopic diffusivity and the macroscopic drift coefficient. This is confirmed in the 
next section by a renormalization group calculation. 
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8 Renormalization Group Equations 



The breakdown of the Ward identity and the absence of any simple relation between the 
effective diffusivity and drift tensors makes it interesting to examine the consequences of 
the renormalization group approach to computing the macroscopic parameters. These 
equations for the diffusivity tensor were written down in refs @, [|]. We complete the 
scheme here with the corresponding equation for the drift coefficient. For simplicity 
we write down the equations for the wavenumber slicing scheme. They are 

dft>ij f d C[ fo(^q A)-D((^) {.9m^"fniQn^nj Qm^TnjQn^ni^Qr^rsQs Qm^miQn^njQr^'rsQs 



dA J (27r) 3 " (q r Krsq s ) 2 

(43) 

and 

dXij f d q . A\rV \ C lmXmiC[nXnjQrXrsQ. 



5(q A)D(q) mi ^ n "-njir '>rs^s (44) 
dA J (2tc) 3 {q r K rs qs) 2 

It is clear from these equations that, as mentioned in a previous paper ||, those 
solutions satisfying boundary conditions for which A°- oc maintain the the propor- 
tionality Ajj(A) oc /%(A) for all A with a constant ratio. 

It is easily checked that of course the isotropic solutions are those of earlier papers, 
namely 

Ky(A) = 4(A) = rcoexp |-I^|A A (0)| 5 tJ (45) 

and 



Atf(A) = Ag(A) = A exp |-i^|A A (0)| ^ (46) 
where 

^=Lwf D ^ • (47 > 

We can examine solutions nearby these isotropic solutions, that are perturbed by a 
small anisotropic change in the local drift tensor. It is no restriction to make this small 
change traceless. We have then 

«*i(A) = 4( A ) + 7 to ( 48 ) 
MA) = Ag(A) + /iy 



where 



^■(oo) = (49) 
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with = . A perturbative analysis of eqs( [l"5| , pHD yields 



2 



If we introduce a variable < s < 1 such that 

d 3 q 



D{q) (51) 



W (2tt) 3 " 

and impose the allowed constraint that r\ = rja = and /i = = we find 



c% 4 Ag 4 A 



(is 15 Kq 15 k 

d/xij 12 A| _ _4_A| 

15^ 15k§ % ' 

These equations are easily integrated and yield the result 



Vij + T7—Vij (52) 



3 k§ J 13 15 k 6XP 1 15 K 2 Q 



Kij = /t exp <^ -- — \ Sij - — — exp <^ - — — \ (53) 



f 1A oU A 4 A o\ f 8 A ol 
A^ = A exp^---2 \6 ij + 1-— -2 ex H~77 71 ^ 



It follows that 



(AK- 1 k = ^f% + ^- exp M^m) • (54) 



K>0 \ J Ao [ 5 

For this near isotropic case at least we see that the proportionality of the macroscopic 
drift and diffusivity tensors is restored exponentially as the renormalization process 
takes place. 

9 Numerical Simulation in the Anisotropic Case 

We tested the above results from the renormalization group calculation against a sim- 
ulation for which the asymmetric local drift tensor is diagonal in the coordinate basis 
with elements 

-a/2 
[i° = | -a/2 | . (55) 
a 
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Ao 


a 




K± 


All 


\± 


1. 


.2 


.701(3) 


.735(3) 


.798(5) 


.660(5) 






(.685) 


(.732) 


(.803) 


(.674) 


1. 


.1 


.698(3) 


.727(2) 


.770(6) 


.687(5) 






(.701) 


(.724) 


(.760) 


(.695) 


1. 


-.1 


.730(5) 


.712(2) 


.672(5) 


.730(5) 






(.732) 


(.709) 


(.673) 


(.737) 


1. 


-.2 


.764(4) 


.701(3) 


.622(6) 


.772(7) 






(.748) 


(.701) 


(.630) 


(.760) 



Table 2: Number of modes in random field is 256. Theoretical values are shown in 
brackets. 



Because of the axial symmetry of the local drift tensor the same property ensures that 
the effective diffusion and drift tensors will have the form 



Ke = «l and A e = \± (56) 





The results of the numerical simulation for certain values of Ao and various values of 
a are shown in Table 2. 

The results with some small discrepancies for the two cases with \oc\ — 2 are in good 
agreement with the predictions of the the renormalization group calculations expressed 
in eqs(^4|) . It is not unreasonable that an asymmetry parameter as large as 20% is the 
limit of applicability of the simple perturbation approach in the previous section. To 
check the results by simulation in finer detail for smaller values of |a] requires higher 
statistical accuracy than can easily be achieved. Our simulations typically involved 256 
particles in 256 velocity fields and required 100/200 processor hours. Nevertheless we 
can conclude with reasonable confidence that the RGC has produces an accurate result 
even in the asymmetric case. It may be that an important condition for the success of 
the RGC is the isotropy of the random field statistics. 



10 Conclusions 

In this paper we have confirmed the importance of the role of drift in the set of effective 
parameters that govern the long time behaviour of a diffusion process in which a particle 
moves subject to molecular diffusion and to the influence of the gradient of a random 
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scalar field. Another way of expressing this is to say that in addition to the the diffusion 
process itself, it is important to consider the effect on the system of long range external 
fields and the strength with which they are coupled to the system. We have exploited 
this effect by computing and measuring the mean drift in a simulation induced by an 
external field of constant gradient. The effects however must also be relevant to all 
external fields of long range. The conclusion, arrived at in previous work 0, H, 3 , 
and confirmed in this simulation is that the effective drift and diffusivity parameters 
are renormalized relative to their local values by identical factors in the isotropic case. 
We showed how the renormalization of the drift parameter could be interpreted as 
a result of biased spatial averaging effects due to the density distribution adopted 
by particles passing through the medium represented by the random field under the 
influence of an external field of constant gradient. 

We also examined a situation in which the isotropy is broken by giving the drift tensor 
an axi-symmetric form. We confirm that the Ward identity suggested in previous work 
as an explanation for the equality of the diffusion and drift renormalization factors 
is indeed broken in lowest order perturbation theory. We do not expect a simple 
relationship between effective diffusion and drift in this case. This is confirmed at 
a theoretical level by using the renormalization group approach in a near symmetric 
situation. The predictions of the theoretical calculation are quite well verified by the 
results of a simulation. The implication is that the renormalization group calculation 
will work reasonably well when the statistics of the random scalar field are isotropic. 
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Figure 1: Vertex diagram 



Figure 2: Drift insertion vertex 



Figure 3: One-loop contribution to E 




Figure 4: One-loop correction to drift insertion 
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Figure 5: One loop correction to vertex 



